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Abstract 

Background: The unfolded protein response (UPR) controls the protein folding capacity of the endoplasmic 
reticulum (ER). Central to this signaling pathway is the ER-resident bifunctional transmembrane kinase/ 
endoribonuclease Irel. The endoribonuclease (RNase) domain of Irel initiates a non-conventional mRNA splicing 
reaction, leading to the production of a transcription factor that controls UPR target genes. The mRNA splicing 
reaction is an obligatory step of Ire 1 signaling, yet its mechanism has remained poorly understood due to the 
absence of substrate-bound crystal structures of Irel, the lack of structural similarity between Irel and other 
RNases, and a scarcity of quantitative enzymological data. Here, we experimentally define the active site of Irel 
RNase and quantitatively evaluate the contribution of the key active site residues to catalysis. 

Results: This analysis and two new crystal structures suggest that Irel RNase uses histidine H1061 and tyrosine 
Y1043 as the general acid-general base pair contributing > 7.6 kcal/mol and 1.4 kcal/mol to transition state 
stabilization, respectively, and asparagine N1057 and arginine R1056 for coordination of the scissile phosphate. 
Investigation of the stem-loop recognition revealed that additionally to the stem-loops derived from the classic Irel 
substrates HAC1 and Xbpl mRNA, I re 1 can site-specifically and rapidly cleave anticodon stem-loop (ASL) of 
unmodified tRNA phe , extending known substrate specificity of Irel RNase. 

Conclusions: Our data define the catalytic center of Irel RNase and suggest a mechanism of RNA cleavage: each RNase 
monomer apparently contains a separate catalytic apparatus for RNA cleavage, whereas two RNase subunits contribute 
to RNA stem-loop docking. Conservation of the key residues among Irel homologues suggests that the mechanism 
elucidated here for yeast Irel applies to Irel in metazoan cells, and to the only known Irel homologue RNase L. 



Background 

The unfolded protein response (UPR) is an intracellular 
signaling pathway that provides homeostatic feedback 
regulation between the endoplasmic reticulum (ER) and 
the gene expression program in the nucleus. To this end, 
the UPR senses the conditions inside the ER, detecting 
an imbalance between newly made proteins and the pro- 
tein folding capacity in the ER, and activates a corrective 
response. For signaling, the UPR uses a transmembrane 
sensor of ER-lumenal unfolded proteins, Irel. 

Irel is an ER membrane-resident receptor that serves 
as a primary signal transduction device in the UPR 
conserved from yeast to mammalian cells [1-4]. 
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Oligomerization of Irel-lumenal domains is thought to 
be a key event in initiating signal propagation across 
the ER membrane that enables the cooperative assem- 
bly of Irel's cytosolic kinase and RNase modules into 
an ordered oligomer with a defined three-dimensional 
structure [5]. The oligomer is stabilized by phosphates 
resulting from autophosphorylation of the kinase 
domain and allows juxtaposition of Irel's RNase 
domains, which presumably activate the RNase [5,6]. 
Once activated, Irel initiates the non-conventional 
splicing of HAC1 mRNA (yeast) or XBP1 mRNA 
(metazoan) by cleaving the mRNA at two conserved 
sites to excise an intron [5,7]. An RNA ligase (tRNA 
ligase in yeast and a still unknown enzyme in 
metazoan cells) rejoins the severed exons to complete 
the reaction. Intron removal allows for the production 



O© 201 1 Korennykh et al; licensee BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative 
BiolVlGCl C^ntrBl Commons Attribution License (http://creativecommons.Org/licenses/by/2.0), which permits unrestricted use, distribution, and 
reproduction in any medium, provided the original work is properly cited. 



Korennykh et al BMC Biology 201 1, 9:47 
http://www.biomedcentral.eom/1 741 -7007/9/47 



Page 2 of 1 5 



of the UPR transcription activators Hacl and XBP1, 
respectively, which upregulate UPR target genes. 

Irel's RNase domain thereby serves the primary role 
in signal transmission. The mechanism underlying acti- 
vation of Irel RNase by oligomerization and the 
mechanism of mRNA recognition and cleavage have 
remained elusive. To date, there are no known struc- 
tural homologues of the Irel RNase domain that could 
help answering these questions. The best attempt at 
defining the mechanism of RNA cleavage has been 
made based on sequence conservation arguments and 
on a crystal structure of Irel dimer with a ligand-free 
RNase domain, which also lacked a fragment of the 
RNase active site [6]. In this work, the RNase dimer was 
proposed to contain two independent catalytic centers, 
one per Irel monomer, which would simultaneously 
accommodate the two RNA stem-loops conserved in all 
known mRNA substrates of Irel [6]. According to the 
proposed model, two stem-loops would form a kissing 
interaction for docking and recognition by Irel. Subse- 
quent studies proposed [5] that such a kissing interac- 
tion is not likely because Irel exhibits no preference for 
RNA substrates with dual stem-loops over substrates 
with a single stem-loop, indicating that a single stem- 
loop structure is the cognate folded substrate of Irel 
RNase. 

It has been also suggested [6] that putative active site 
residues are positioned similarly in Irel RNase and in 
pre-tRNA splicing endonuclease (SEN), although these 
two endoribonucleases share no apparent sequence or 
structural homology and cleave dissimilar RNA sub- 
strates (stem-loop versus bulge-helix-bulge motif). The 
authors point out [6] that the putative catalytic residues 
in Irel do not align completely with those in SEN and 
would require a 5-A displacement of Y1043 (in yeast 
Irel) to bring it into a position corresponding to that of 
Y249 in SEN. A 5-A conformational change in Irel 
RNase upon substrate binding has been proposed based 
on these arguments [6]. A new cocrystal structure of 
Irel RNase with an oligonucleotide bound, which we 
provide in this work, does not support considerable con- 
formational changes in the position of Y1043. Therefore, 
Irel and SEN apparently cleave RNA using different 
arrangements of active site residues. Underscoring this 
difference, substitution of a single catalytic histidine 
residue H1061 produced orders-of-magnitude greater 
effect on the catalytic activity of Irel (this work) com- 
pared to only a 28-fold rate reduction in SEN [8]. 

Conservation of residue H1061, the only invariant his- 
tidine in the RNase domain of Irel, was used to propose 
that H1061 marks the catalytic center and serves for 
general acid-general base catalysis [6], In agreement 
with this model, mutation H1061A reduced RNA clea- 
vage [6]. However, this mutation was analyzed using a 



qualitative approach that had a narrow dynamic range 
and did not distinguish between effects on binding and 
catalysis. Furthermore, mutation H1061A not only 
removed the ability of histidine to transfer protons but 
also the ability to form hydrogen bonds, which alone 
could explain the observed effect on RNA cleavage, 
without invoking H1061 in general acid-general base 
catalysis. Considering that ribonucleases can cleave RNA 
stem-loops by using for proton transfer basic amino 
acids instead of conventional ionizable residues (histi- 
dine, glutamate and aspartate) [9], the very involvement 
of H1061 in general acid-general base catalysis could 
not be deduced convincingly from any of the available 
data. 

To understand the mechanism of RNA cleavage by 
Irel, we combined quantitative analysis of rationally 
designed Irel mutants and X-ray crystallography. We 
show that proton transfer by histidine H1061 contri- 
butes greater than five orders of magnitude (> 7.6 kcal/ 
mol) to catalysis of RNA cleavage, which experimentally 
defines a key catalytic functionality in the active center 
of Irel RNase. Using the non-disruptive Y1043F muta- 
tion we experimentally demonstrate a catalytic role of 
the OH group of Y1043 and also show an important 
role in RNA binding of a helix-loop element (HLE) resi- 
due, R1039. Our findings suggest parallels between cata- 
lytic centers of Irel and the well characterized RNases 
A and Tl and propose an unexpected mechanism of oli- 
gomerization-induced in trans RNA recognition by Irel. 

Results 

Computational analysis of the RNA binding properties of 
Irel 

To glean insights into Irel interactions with HAC1I 
XBP1 mRNA, we analyzed the RNA binding properties 
of the entire surface of the Irel oligomer [5], Surface 
electrostatics calculations [10] are commonly used for 
such analyses, but they produce ambiguous results 
because electrostatics is not the sole determinant of 
RNA binding. Here we used a modified optimal protein- 
RNA area (OPRA) algorithm that previously demon- 
strated high prediction accuracy [11]. We refined this 
method by adding phylogenetic conservation to the 
scores and termed the resulting algorithm 'conservation- 
weighted OPRA' (cwOPRA; to be described elsewhere). 
The cwOPRA analysis of RNA binding properties of the 
Irel oligomer produced a high-contrast image of the 
surface RNA binding propensity. A single RNA binding 
site was strongly predicted in the center of every back- 
to-back RNase dimer [6]. The RNA binding site is prox- 
imal to the dynamic helix-loop element HLE, which 
becomes more ordered upon formation of Irel/Irel 
interfaces in the high-order oligomer and has been pro- 
posed to serve as a part of the active site [5] (Figure la, 
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Figure 1 Oligonucleotide binding to Irel RNase (a) Ire 1 oligomer (PDB ID 3fbv) colored by domains (top) and RNA-binding propensity 
(bottom), (b) Close view of symmetric back-to-back RNase dimer from the oligomer in (a) with electron density for dCdCdGdCdAdG from the 
C222 crystal structure superimposed. Electron density colocalizes with the area of the strongest RNA binding propensity and is in direct contact 
with conserved histidine H1061 and the helix-loop element (HLE). (c) Binding of an RNA stem-loop HP21 and of dCdCdGdCdAdG to Irel KR32 
measured via inhibition of 32 P-5'-HP21 cleavage by unlabeled HP21 (K m = 5 ± 1 uM) or dCdCdGdCdAdG (K; app = 200 ± 50 uM). Reactions 
contained 20 mM 4-(2-hydroxyethyl)-1-piperazine-ethanesulfonic acid (HEPES) (pH 7.4), 70 mM NaCI, 5% glycerol, 2 mM Mg(CH 3 COO) 2 , 2 mM 
ADP, 4 mM dithiothreitol, 3 uM total lre1KR32 and < 1 pM 32 P-5'-HP21, and were conducted at 30°C. Measurements of HP21 binding at 1 uM 
and 10 uM of total Irel KR32 are shown in Additional file 1, Figure SI. (d), Wall-eyed stereoview of non-crystallographic symmetry (NCS)-averaged 
F G -F C electron density map at the RNase active site. Contour levels are 7a (dark blue) and 13c (red). Side chain positions are from the crystal 
structure 3fbv. 
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b). The rest of the Irel surface was predicted not to 
interact with RNA (we note, however, that the crystal 
structure is lacking a highly basic N-terminus, which is 
disordered and could serve as a secondary RNA binding 
site). We next validated this computational prediction 
using a cocrystal structure of Irel with an 
oligonucleotide. 

Crystal structure of Irel with dCdCdGdCdAdG bound to 
the RNase domain 

For crystallization trials, we tested a panel of RNase- 
resistant RNA analogs with 2'-deoxy or 2'-methoxy sub- 
stitutions throughout and a variety of Irel constructs. 
Only a single combination led to formation of diffract- 
ing crystals as discussed below. Attempts to soak mono- 
nucleotides, dinucleotides, trinucleotides or 
oligonucleotides (0.5-20 mM) into previously character- 
ized Irel crystals [5] were unsuccessful. No additional 
difference density was observed, which could correspond 
to the soaked-in substrate analogs. 

By contrast, cocrystallization of IrelKR32 with a 6- 
mer 2'-deoxy-substituted oligonucleotide bearing the 
consensus sequence of an Irel splice site CCGCAG 
(the use of 2'-deoxy substitutions prevents degradation 
by Irel RNase) resulted in crystals of the 
IrelKR32-dCdCdGdCdAdG complex. IrelKR32 con- 
struct contains the Irel's kinase domain, the RNase 
domain, and a 32-amino-acid linker previously shown 
to be required for optimal RNase activity [5]. The 
dCdCdGdCdAdG oligonucleotide bound to IrelKR32 
with Kj app = 200 uM determined from inhibition of the 
IrelKR32-catalyzed cleavage of HP21, a small stem/ 
loop structure containing the 3' splice site of XBP1 
mRNA [5]. This value is approximately 40-fold weaker 
than that for the cognate stem-loop (K m = 5 uM, Fig- 
ure lc and Additional file 1, Figure SI) possibly due to 
a more flexible structure of the single-stranded oligo- 
nucleotide dCdCdGdCdAdG compared to HP21 RNA 
stem-loop and/or due to the lack of 2'OH groups. 
Unlike the oligonucleotide-free Irel, which crystallizes 
in space group P2!2!2 [5], the IrelKR32 dCdCdGdC- 
dAdG oligomer crystallized in space group C222 and 
diffracted to 6.6 A (Additional file 1, Table SI). The 
presence of the oligonucleotide allowed IrelKR32 to 
crystallize without a ligand occupying the kinase ATP 
binding pocket, providing an important insight into 
the mechanism of cofactor-dependent activation of 
Irel [12]. 

We used the 3.2-A structure of Irel crystallized in the 
absence of a nucleic acid substrate [5] for molecular 
replacement and subsequent rigid-body refinement, as 
described in Methods and in the accompanying manu- 
script [12]. At 6-7 A resolution, the secondary structure 
of proteins and the sugar-phosphate backbone, which is 



the most electron-dense part of nucleic acids, can be 
visualized in ribonucleoprotein complexes [13,14]. 
Accordingly, we observed a strong F 0 -F c difference elec- 
tron density at the core of the RNase domain, which 
coincided with the cwOPRA-predicted RNA binding site 
(Figure lb). The electron density was delineated by two 
separate contours related by twofold symmetry (Figure 
lb, d). Each contour (at level 7a) defines an elongated 
volume of approximately 20 A in length, consistent with 
the length of one dCdCdGdCdAdG molecule bound per 
RNase monomer. We placed Irel side chains according 
to the 3.2-A-resolution search model used for the mole- 
cular replacement (Figure Id). Such an approximation is 
valid because the protein backbone traces and the qua- 
ternary structure of Irel oligomers superimpose in the 
6.6A and the 3.2A structures, even though they crystal- 
lized in different space groups (Pl^^ vs C222). Based 
on the complete superposition of the polypeptide back- 
bones in IrelKR32APY29 [5] and in Irel-dCdCdGdC- 
dAdG complexes, we can rule out the possibility of a 
large-scale movement in the catalytic core of Irel upon 
binding of the substrate analog. 

Probing the catalytic mechanism of Irel RNase via 
mutagenesis: role of H1061 

The electron density for dCdCdGdCdAdG is located in 
the cleft formed by the HLEs [5]. The position of the 
electron density suggests that eight protonatable resi- 
dues of Irel located proximal to the substrate mimic 
could, in principle, contribute to general acid-general 
base catalysis leading to the formation of the 2', 3'-cyclic 
phosphate [15]. Because Irel does not require divalent 
metal ions for RNase activity [5], it must rely exclusively 
on amino acid side chains to carry out catalysis. In the 
identified group of amino acids, only H1061 and D1064 
have pK a values near neutrality, analogous to the cataly- 
tic residues involved in proton transfer in well charac- 
terized ribonucleases, such as RNase A [16] and RNase 
Tl [17]. The H1061 is invariant in all known Irel 
sequences and the only Irel homologue RNase L [18]. 
By contrast, D1064 is mutated in several organisms to 
residues incapable of proton transfer (Additional file 1, 
Table S2), indicating that it is dispensable for catalysis. 

Mutation of the catalytic histidine to alanine, aspara- 
gine or glutamine diminishes catalysis by ribonucleases 
Tl and A by 10 2 to 10 4 -fold [17,19,20]. In a previous 
study, an H1061A mutation in Irel was reported to 
inhibit the RNase activity; however, the quantitative 
effect of this mutation on catalysis has not been mea- 
sured [6]. To determine the quantitative contribution of 
H1061 to catalysis, we expressed and purified a mutant 
of IrelKR32 in which H1061 was replaced by aspara- 
gine. We chose mutation to asparagine because it best 
preserves the space filling and hydrogen bonding 
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properties of histidine, while removing histidine's ability 
to transfer protons and to support catalysis. We 
observed that IrelKR32(H1061N) displays a > 3 x 10 s - 
fold rate reduction compared to wild-type IrelKR32 
(Figure 2a), strongly indicating that H1061 plays a cen- 
tral part in catalysis. 

Next, we ruled out a possibility that the large effect of 
the H1061N mutation arose from a perturbation of oli- 
gomerization. To this end, we monitored the optical 
density of IrelKR32(H1061N) solutions at 500 nm, 
which provides a quantitative measure of IrelKR32 oli- 
gomerization (Methods and [12]). IrelKR32(H1061N) 
mutant exhibited an indistinguishable oligomerization 
profile from that of wild-type IrelKR32 (Figure 2b), sug- 
gesting that the H1061N mutation does not interfere 
with oligomer formation. 

To further ascertain that the dramatic rate reduction 
does not arise from more local conformational perturba- 
tions of Irel RNase, uncoupled from oligomerization, we 
crystallized IrelKR32(H1061N) in complex with the syn- 
thetic ligand APY29, which binds in Irel's kinase ATP 
binding pocket [5]. We determined a 3.65-A-resolution 
crystal structure of this mutant (Additional file 1, Table 
S3), which proved identical to that of wild-type 
IrelKR32 (Figure 2c; Additional file 1, Figure S2). Our 
functional and structural results pertaining to H1061N 
mutant therefore indicate that the greater than 3 x 10 - 
fold catalytic impairment of IrelKR32(H1061N) arises 
from a direct role of H1061 in proton transfer during 
RNA cleavage. 

Due to the twofold symmetry of the Irel RNase dimer 
[6], a single H1061N mutation removes two H1061 resi- 
dues in the vicinity of the RNase active site (Figure lb). 
We thus wondered whether the catalytic impairment of 
the H1061N mutant might arise from this symmetry: if 
Irel required both H1061 residues together to catalyze 
phosphodiester cleavage, we would de facto have mea- 
sured activity of a double mutant. To explore this possi- 
bility, we built Irel oligomers in which each catalytic 
H1061 would be paired with a non-catalytic H1061N in 
chimeric RNase dimers. To this end, we titrated catalyti- 
cally inactive IrelKR32(H1061N) into a solution of wild- 
type IrelKR32 held at a constant concentration just 
below its oligomerization threshold. For enzymes that 
follow Michaelis-Menten kinetics, addition of a catalyti- 
cally inactive mutant to a wild-type reaction inhibits the 
rate, as inactive enzyme sequesters the available sub- 
strate. By contrast, we observed approximately tenfold 
activation of IrelKR32 before onset of inhibition at 
excessive concentrations of IrelKR32(H1061N) (Figure 
2d). Stimulation occurs because wild-type IrelKR32 
becomes trapped in a wild-type/H1061N hetero-oligo- 
mer as enzyme concentration increases and causes Irel 
to oligomerize. Kinetic modeling of this experiment 



shows that such a response is expected only if wild-type: 
H1061N heterodimers in the IrelKR32 oligomer were 
active RNases (Additional file 1, Figure S3 and Addi- 
tional file 2). We conclude that only one of the two 
H1061 residues in the twofold symmetric RNase dimer 
is required for RNA cleavage. 

Proximity to the electron density of the substrate, 
mutagenesis analysis, evolutionary conservation and the 
neutral pK a value of histidine all converge on the sug- 
gestion that H1061 serves a catalytic function. Cleavage 
of an RNA backbone by other non-metal RNases 
involves pairs of conserved residues that carry out gen- 
eral acid-general base catalysis and lie at a fixed mutual 
separation of 6-8 A (for example, PDB ID lr5c, lrga; 
Figure 3a). In Irel, only two invariant amino acids, 
R1056 and Y1043, are located 5.5 A and 8 A from 
H1061, respectively (Additional file 1, Table S2), and 
hence emerge as candidates to be a catalytic partner of 
H1061. Steric considerations discussed below and high 
pK a value of arginine argue that R1056 is unlikely to ful- 
fill this role, suggesting Y1043 as the most plausible 
partner of H1061 (Figure 3b, c; see Discussion). We 
note that the R1056 side chain is deeply buried in the 
active site, which would make it difficult to interpret 
results of mutational analyses. We therefore did not 
pursue analysis of R1056 at this time. 

Probing the role of Y1043 as the catalytic partner of 
H1061 

Previously, it was suggested that if Irel RNase resembled 
SEN endoribonuclease, then tyrosine Y1043 may serve 
as the general base [6]. The fact that Y1043 resides 
approximately 5 A from the analogous tyrosine in SEN 
was used to propose that Irel RNase should undergo an 
obligate conformational change upon binding of sub- 
strate to position Y1043 relative to H1061 as in SEN. 
To test this possibility, we examined the electron density 
for the Y1043-bearing loop in the IrelKR32 dCdCdGdC- 
dAdG complex (Figure 3b; Additional file 1, Figure S4). 
We found that backbones of the HLE and of the adja- 
cent Y1043-bearing loop occupy the same position as in 
the 3.2-A crystal structure obtained without bound oli- 
gonucleotide (PDB ID 3fbv). These findings support a 
model that tyrosine Y1043 does not move over a signifi- 
cant distance upon substrate binding. 

To assess whether Y1043 contributes to catalysis, we 
prepared and quantitatively characterized the mutant 
IrelKR32(Y1043F), which is predicted to disable any 
involvement of Y1043 in proton transfer while ensuring 
minimal structural perturbation to the protein compared 
to a potentially far more disruptive Y1043A mutation 
tested previously [6]. We obtained an activation profile 
with IrelKR32(Y1043F) (Figure 3d, open circles). At 
enzyme concentrations exceeding 3 uM, the observed 
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Figure 2 Probing the role of H1061 in catalysis, (a) Effect of H 1 061 N mutation on the RNase activity of Irel. Error bars show standard error of 
single-exponential fitting, (b) Oligomerization profiles of Ire1 KR32 and three mutants, Irel KR32(H1061 N), Irel KR32(Y1043F) and lre1KR32(R1039A) 
measured using light absorbance at 500 nM (see Methods). Measurements were conducted at room temperature using the same reaction 
buffers as in Figure 1c. (c) Comparison of Irel structures from oligomers formed by wild-type Irel (PDB ID 3fbv, space group P2 1 2 1 2) and by 
H 1 061 N mutant (new 3.65A crystal structure, also in space group P2,2 1 2). Shown are overall oligomer architectures (left), 2F 0 -F C electron density 
map for the H1061N RNase domain contoured at 13a (middle) and superposition of the RNase domains (right). Simulated annealing F D -F C omit 
map for H1061N calculated without non-crystallographic symmetry (NCS) is shown in Additional file 1, Figure S2. (d) In trans activation of 
re1KR32 by the non-catalytic mutant H1061N. Concentration of wild-type Irel KR32 was 1.5 uM, concentration of the H1061N Irel varied 
between 0.1-15 uM. Arrow marks the point of equivalent concentrations. Reactions were conducted under single-turnover conditions using the 
same reaction buffer as in Figure 1c. 
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Figure 3 Catalytic mechanism of Ire1 RNase. (a) Catalytic residues of RNase T1 and RNase A with productively bound nucleotides (PDB ID 
1r5c and 1rga). (b) Left panel: simulated annealing (1000 K) omit-electron density map F Q -F C for the helix-loop element (HLE; green) and 
adjacent loop containing Y1043 and the helix a4 of RNase domain (contoured at 5c). Refinement was conducted with residues 1032-1057 
deleted from all 7 monomers in the asymmetric unit, without using non-crystallographic symmetry (NCS). An alternative view of the electron 
density (after rotation around vertical axis) is shown in Additional file 1, Figure S4. Right panel: position of the scissile phosphate in the active 
site of Ire1 RNase. The strongest peak of the electron density (F 0 -F 0 contoured at 12c) was used to position the phosphate. Side chains were 
taken from PDB ID 3fbv). (c) Two plausible catalytic mechanisms in the active site of Ire1 RNase. (d) Enzyme titration profiles for lre1KR32, 
Irel KR32(Y1043F) and for lre1KR32(R1039A) under single-turnover conditions. Reactions were conducted as in Figure 1c. 
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rate constant for cleavage of HP21 RNA plateaued at a 
first-order single-turnover rate constant k 2 = 0.022 s' . 
This value is approximately tenfold smaller than the 
corresponding value for wild-type IrelKR32 (k 2 = 0.25 s" 
1 ), demonstrating that the OH group of Y1043 is impor- 
tant for catalysis. Similar to the H1061N mutation, the 
Y1043F mutation had no effect on the oligomerization 
equilibrium of IrelKR32 (Figure 2b). The deleterious 
effect of the Y1043F mutation was greater at subsaturat- 
ing k 2 /Ki/ 2 conditions encountered below approximately 
3 uM enzyme concentrations, indicating that Y1043 also 
contributes to binding of the RNA substrate (Figure 3d). 

Probing the role of the HLE in RNA recognition 

Location of the dCdCdGdCdAdG electron density in 
proximity to the HLE (Figure lb, d) confirms the pre- 
viously posed hypothesis [5,6] that this element serves as 
an important part of the active site. It was proposed that 
residues from HLEs could serve for binding and recogni- 
tion of RNA substrates, rather than for catalysis. Muta- 
tions in the HLE are therefore expected to provide a 
valuable control to contrast the catalytically important 
mutations of H1061 and Y1043 described above. To test 
the role of the HLE in RNA cleavage, we prepared 
IrelKR32 with a mutation of an arginine residue within 
the HLE, IrelKR32(R1039A) and measured the titration 
rate profile for the IrelKR32(R1039A) mutant (Figure 3d, 
diamonds). At saturating enzyme concentrations 
IrelKR32(R1039A) achieved the same first-order single- 
turnover rate constant k 2 as did wild-type IrelKR32 (Fig- 
ure 3d, closed circles), indicating that R1039 makes no 
contribution to stabilization of the transition state. Under 
subsaturating conditions, however, the R1039A mutant 
became deleterious, revealing a role of R1039 in substrate 
binding. Analogous to other mutations tested here, 
IrelKR32(R1039A) exhibited an unperturbed oligomeri- 
zation profile (Figure 2b), indicating that this mutation 
too does not affect the enzyme's oligomerization equili- 
brium. Both structural (Figure Id) and functional data 
with IrelKR32(R1039A) therefore show that the HLE 
plays an important role in recognition of RNA substrates. 

Electron density-guided docking of an RNA stem-loop 
into Irel RNase (below) suggests that the HLE of one 
RNase monomer may contribute to RNA cleavage in the 
active site of an adjacent RNase monomer, acting in 
trans. To test this possibility, we performed in trans 
activation experiments as in Figure 2d, in which we ana- 
lyzed the activity of IrelKR32(wild-type)/IrelKR32 
(R1039A) heterodimers. To ensure that effect of the 
R1039A mutation is indeed measured in trans, we addi- 
tionally disabled the catalytic site of the R1039A mutant 
by the H1061N mutation, as outlined in Figure 4a. The 
experiment was conducted as in Figure 2d, but con- 
tained lower concentration of wild-type IrelKR32 (0.3 



uM or approximately 0.06 K m ) to ascertain a regime in 
which subsaturating conditions (k 2 /Ki/2) were main- 
tained during the titration. This was an important con- 
sideration because if the reactions were allowed to 
approach the first-order k 2 regime, the deleterious effect 
of the R1039A mutation on binding would have been 
masked (note loss of effect from the R1039A mutation 
in Figure 3d when the reaction plateaued at E 0 > 5 uM). 
Using this experiment, we found that the RNase activity 
of wild-type IrelKR32 stimulated in trans by the double 
mutant IrelKR32(R1039A, H1061N) was an order-of- 
magnitude lower than with the H1061N single mutant 
(Figure 4b). The chimeric dimer with the double muta- 
tion (G in Figure 4a) therefore binds RNA more weakly 
than the chimeric dimer with the single mutation (E in 
Figure 4a), whereas the only difference between these 
two cases is an intact HLE (E) versus a HLE bearing a 
single R1039A mutation (G). It follows that the RNA 
cleavage reaction catalyzed by one RNase monomer is 
sensitive to a mutation in the HLE of the paired mono- 
mer. These observations strongly support a model of a 
composite RNA-binding surface in Irel RNase, which 
also emerges from the electron density-guided RNA 
stem-loop docking described below. 

Irel site-specifically cleaves stem-loop anticodon of 
tRNA phe 

Irel recognizes and site-specifically cleaves conserved 
sites in HAC1 and XBP1 mRNA, which are character- 
ized by a stem containing approximately five Watson- 
Crick base pairs and a seven-residue loop with a con- 
sensus sequence CNGNNGN [21]. To understand how 
Irel might recognize such stem-loops, we searched 
bioinformatically for RNA heptaloops with similar fea- 
tures. This search showed that the anticodon stem-loop 
of tRNA phe closely resembles the Irel consensus sub- 
strate (Figure 5a). To verify that this sequence similarity 
meaningfully identifies an Irel substrate, we conducted 
RNA cleavage assays and compared the kinetics of clea- 
vage by Irel of tRNA phe anticodon stem-loop and clea- 
vage of the stem-loop derived from XBP1 mRNA, HP21. 
Irel cleaved both stem-loops in a single site and with 
similar rates (Figure 5b), indicating that the tRNA phe 
anticodon stem-loop indeed provides a good model of 
the cognate Irel stem-loops from XBP1 and HAC1 
mRNA. As multiple crystal structures of tRNA Phe have 
been already solved, this finding enabled us to model 
the docking of an RNA stem-loop to the Irel RNase 
active site (Figure 5c-e). 

Discussion 

The architecture of Irel RNase 

Our conclusion that each Irel monomer contains a 
complete catalytic center in the context of the Irel 
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Figure 4 /n trans activation of wild-type Ire1 RNase by Ire1 with RNase with a double mutation (a) Schematic representation of RNase 
dinners formed via interface IF1 C [5], Due to a twofold symmetry, each RNase dimer can accommodate RNA stem-loop in two equivalent 
orientations (such as A and B). Both orientations are productive for wild-type Irel. Neither orientation is productive for the H1061N mutant (C) 
and (D) due to disrupted catalysis. Only one orientation (E) is productive for the wild-type/H1061N chimera. Only one orientation (G) is also 
productive for the wild-type/H1061N+R1039A chimera, however this orientation is impaired by the R1039A mutation in the helix-loop element 
(HLE). (b) Titration of wild-type Irel KR32 with rel KR32(H 1 061 N) single mutant and Ire1 KR32(H1 061 N, R1039A) double mutant. Reactions were 
conducted as in Figure 2d but contained 0.3 uM wild-type lre1KR32 throughout the titrations to ensure subsaturating (k 2 /K 1/2 ) regime for RNA 
cleavage. 



oligomer agrees with the model proposed previously 
based on experiments in which several Irel variants 
with dissimilar RNase-inactivating mutations were 
mixed [6]. The previous observation showed that none 
of the mutant mixtures have RNase activity and was 
interpreted in favor of individual active sites over a com- 
posite catalytic center (created by association of two 
RNase monomers) [6]. Importantly, however, this lack 



of RNase activity can be reasonably reconciled with 
both, per monomer and composite active site architec- 
tures, as well as with a third possibility of indirect 
defects in Irel RNase due to misfolding or conforma- 
tional perturbation of the mutants. By contrast, the 
observation of electron density for a substrate mimic in 
each RNase monomer (Figure lb, d) and in trans activa- 
tion of wild-type IrelKR32 by IrelKR32(H1061N) 



Korennykh et al. BMC Biology 201 1, 9:47 
http://www.biomedcentral.eom/1 741 -7007/9/47 



Page 1 0 of 1 5 



a b 




substrate 



product 



HP21 tPhe 



C 

tRNA phe «MiaA complex 




anticodon 




tRNA Phe 






non-productive 




monomer A 



monomer B 



1. catalysis 



2. recognition 




HLE (B) 



stem-loop 



Figure 5 Model of stem-loop recognition via a composite protein surface (a) Stem-loops derived from XBP1 mRNA (HP21) and 
phenylalanine tRNA phe (tPhe). Loop sequence differences are colored yellow in tPhe. (b) Time courses (0-5 min) for cleavage of 5'- 32 P-labeled 
stem-loops HP21 and tPhe by Irel KR32 (3 uM) analyzed by polyacrylamide gel electrophoresis, (c) Superposition of crystal structures of tRNA phe 
(PDB ID 1ehz) and tRNA phe from tRNA phe -MiaA complex (PDB ID 2zm5). Black arrow shows direction of nucleophile attack, orange arrow shows 
direction of leaving group, (d) Manual rigid-body docking of the base-flipped tRNA phe anticodon into Irel crystal structure. Electron density for 
the scissile phosphate (contoured at 12c) and the proposed catalytic mechanism (Figure 3c-l) were used as docking constraints. P designates 
scissile phosphate. The helix-loop element (HLE) is colored green, (e) Surface rendering of the model in (d). The model suggests cleavage of a 
stem-loop by one active site and recognition by a composite RNA binding surface. 
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(Figure 2d) provide strong experimental support for an 
individual catalytic center in each RNase monomer. 

Position of the scissile phosphate and the catalytic 
mechanism of Irel RNase 

We deduce the position of the scissile phosphate of a 
bona fide Irel RNA substrate in this catalytic center 
using the structure of the IrelKR32 complexed with a 
DNA hexanucleotide, dCdCdGdCdAdG, which mimics 
the cognate Irel substrate. This substrate analog was 
chosen due to present experimental limitations imposed 
by the inability of Irel to cocrystallize with other oligo- 
nucleotides. While we recognize that the deoxy-substi- 
tuted single-stranded substrate may be structurally 
distinct from RNA stem-loops, we surmise that its bind- 
ing site on Irel represents that of the HP21 RNA 
because (i) dCdCdGdCdAdG efficiently inhibits cleavage 
of HP21 by IrelKR32, and (ii) the electron density in 
the RNase active site directly contacts the critical cataly- 
tic residue H1061, which is involved in proton transfer 
and contributes > 3 x 10 5 -fold to catalysis. At high con- 
tour levels, we could identify three near-spherical peaks 
in a F 0 -F c Fourier map. The central peak was strongest 
and was visible in most Irel monomers in the asym- 
metric unit, while the other two peaks could be 
observed in several but not all Irel monomers. Because 
the central peak most likely represents the most-ordered 
phosphate group of the oligonucleotide and is located 
between the functionally validated catalytic residues 
H1061 and Y1043 (Figures 2a and 3d), we suggest that 
it corresponds to the position of the scissile phosphate 
of the substrate. 

Placement of this phosphate and comparison of the 
active site of Irel with the active sites of RNase Tl and 
A suggest putative roles of Irel residues in catalysis 
(Figure 3a, b). Similar to Irel, ribonucleases Tl and A 
do not use metal ions [5] and produce 2',3'-cyclic phos- 
phate [21]. Importantly, the mechanisms of RNA clea- 
vage have been established in detail for both RNase Tl 
and RNase A [20,22-24]. 

Structural comparison shows that Irel and RNases Tl 
and A have a single positively charged residue contact- 
ing the phosphate: R1056 (Irel), R77 (RNase Tl) and 
K41 (RNase A). These positively charged residues serve 
for catalysis by stabilizing the negative charge buildup 
on non-bridging oxygens in the transition state, as has 
been shown for RNase A [20]. By analogy, R1056 of Irel 
may also participate in transition state stabilization and 
is unlikely to directly participate in proton transfer. The 
last notion is supported by steric considerations: R1056- 
phosphate-H1061 connectivity would form an atypically 
acute angle for a general acid-phosphate-general base 
triad, placing phosphate entirely off the line connecting 
the two protein residues implicated in proton transfer 



(compare Figure 3a and 3b). This arrangement would 
not be compatible with in-line geometry required for 
RNA cleavage to form a 2',3'-cyclic product. 

In RNase Tl and RNase A, a single uncharged hydro- 
philic residue is positioned off the line connecting gen- 
eral acid, scissile phosphate and general base, and at a 
hydrogen-bonding distance from the phosphate oxygens: 
Y38 in RNase Tl and Nil in RNase A, respectively. In 
RNases Tl and A, these residues play a catalytic role by 
stabilizing the transition state via coordination of non- 
bridging oxygen [25]. Position of the electron density in 
Irel (Figure 3b, right panel) predicts an analogous cata- 
lytic function for N1057. 

Finally, all three RNases contain two residues that act 
as a general acid/general base pair. For RNases Tl and 
A, general acids are H92 and HI 19 at a distance of 3.7 
and 3.8 A from the scissile phosphate, respectively. Gen- 
eral bases are H40 and H12 at a distance 5.6 A and 3.8 
A, respectively. The position of the electron density for 
the phosphate between the H1061 and Y1043 (Figure 
3b) agrees with the role of these residues in Irel as the 
general acid/general base pair (Figure 3c). Whereas the 
role of H1061 parallels that of the histidines in RNase 
Tl and A, a catalytic tyrosine is unusual because of the 
high pK a of approximately 10. Experimental testing of 
the catalytic role of Y1043 using an Y1043F mutation 
revealed that the hydroxy group of Y1043 contributes a 
rate enhancement of approximately tenfold (or approxi- 
mately 1.4 kcal/mol) to catalysis. Our data therefore are 
in agreement with the catalytic role of Y1043 and sug- 
gest two kinetically equivalent catalytic mechanisms 
shown in Figure 3c. We note that the tenfold transition- 
state effect of the Y1043F mutation is modest and thus 
it could also be consistent with a third possibility that 
the OH group of Y1043 plays an indirect role in cataly- 
sis and forms an important hydrogen-bonding network 
in the active site rather than participates in proton 
transfer directly. This alternative possibility cannot be 
ruled out until an atomic resolution cocrystal structure 
of an Irel-RNA complex becomes available. If Y1043 
serves for proton transfer directly, then its relatively 
small catalytic contribution compared to the > 3 x 10 - 
fold contribution of the partnering H1061 residue or the 
10 to 10 -fold contributions of the catalytic histidines 
in RNase A or Tl [17,20] could be accounted for quan- 
titatively by the 10 3 to 10 4 higher pK a of tyrosine com- 
pared to histidine. 

While our data support the roles of the H1061 and 
Y1043 as a general acid/general base catalytic pair, they 
do not unambiguously define which residue serves as 
the general acid or the general base. Such an assignment 
has proven challenging for RNases and ribozymes, even 
when high-resolution cocrystal structures with substrate 
analogs or products were available [20,26]. At present, 
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we cannot formally distinguish between the two plausi- 
ble mechanisms for RNA cleavage by Irel shown in Fig- 
ure 3c. As we argue below, the scheme shown in the 
left panel appears more compatible with steric consid- 
erations upon docking an RNA stem-loop in Irel RNase. 

Our model suggests that the Irel RNase active site is 
unique and does not superimpose with other structu- 
rally characterized RNases. However, the arrangement of 
the catalytic residues in Irel and their interactions with 
the scissile phosphate resemble those in the active sites 
of ribonucleases Tl and A. By extension, we predict 
that the same mechanism applies to mammalian Irel 
and RNase L, the Irel homologue involved in the inter- 
feron response pathway, in which all relevant residues 
discussed here for Irel are strictly conserved. 

A possible mechanism of stem-loop RNA recognition by 
Irel RNase 

The crystal structure of tRNA phe (PDB ID lehz) shows 
the stem-loop in a conformation that would preclude 
cleavage by Irel (Figure 5c, yellow trace). The RNA clea- 
vage reaction imposes strict stereochemical requirements 
such that an incoming nucleophile can attack the phos- 
phate in line with the leaving group. A similar stereoche- 
mical problem was encountered in studies of sarcin/ricin 
GAGA tetraloop cleavage by ribonucleases a-sarcin and 
restrictocin. A cocrystal structure of restrictocin with the 
sarcin/ricin stem-loop RNA suggested that scissile bases 
should undergo base flipping to position the incoming 
2'OH in the 2'-endo conformation and to allow for in- 
line attack during catalysis [27] . Guided by this concept, 
we searched the PDB database for a tRNA phe structure, 
in which the loop guanosine at position 3 is flipped out 
such that its ribose ring adopts the 2'-endo conformation. 
We found that the crystal structure of tRNA phe com- 
plexed with the tRNA modifying enzyme MiaA [28] 
serves as a suitable model (Figure 5c, blue trace). 

Accordingly, we used the stem-loop from the tRNA- 
phe -MiaA structure for rigid-body docking into Irel 
active site. As a restraint, we placed the scissile phos- 
phate into the strongest electron density peak (Figure 
5d) and the rest of the loop to allow sterically permissi- 
ble interactions expected from the cleavage mechanism 
(I) in Figure 3c. The docked model has no backbone 
intersections or clashes of backbone atoms. By contrast, 
cleavage mechanism (II) could not be satisfied using this 
docking approach due to seemingly irresolvable inter- 
chain clashes. The resulting model predicts that only 
one stem-loop can bind to the Irel RNase dimer and 
that one RNase domain in the dimer cleaves RNA, 
whereas the other RNase monomer completes the RNA- 
contacting surface in trans, suggesting that both RNase 
monomers and both HLEs participate in recognition of 
the RNA stem-loop structure (Figure 5e). 



Our model that the HLE of one RNase monomer con- 
tributes to RNA cleavage by the catalytic center of a dif- 
ferent RNase monomer is functionally supported by the 
in trans complementation assay described above (Figure 
4). This experiment demonstrates that cleavage by the 
active site of an RNase monomer is impaired upon 
removal of an RNA binding residue from an adjacent, 
non-catalytic monomer. Therefore, although Irel RNase 
has individual per-monomer catalytic centers, two 
RNase monomers are involved in RNA cleavage and 
apparently form a composite RNA binding surface cre- 
ated via dimerization. The high-order oligomerization of 
Irel completes the active site by stabilizing the highly 
dynamic HLE [5] in a conformation suitable for recogni- 
tion of cognate RNA stem-loops. This mechanism could 
explain how self-assembly of Irel into an oligomer acti- 
vates the RNase function. 

Irel cleavage of the tRNA phe anticodon loop 

We found that, strikingly, Irel can recognize and cleave 
the anticodon stem/loop structure of tRNA phe in addi- 
tion to the cognate splice site stem-loops of HAC1 and 
XBP1 mRNA. This observation redefines the known 
substrate specificity of Irel and suggests that additional 
stem-loops besides consensus sites in HAC1 and XBP1 
mRNA could be cleaved by Irel RNase. Whether 
tRNA phe is an in vivo substrate of Irel remains to be 
investigated. In principle, tRNA destruction could serve 
to reduce the cell's translational capacity and thereby 
reduce protein influx into the ER, akin to protein 
kinase-like ER kinase (PERK) activation in metazoan 
cells. However, the scissile guanine in the anticodon of 
tRNA phe carries a post-transcriptional 2'-methoxy modi- 
fication [29], which precludes cleavage of the modified 
tRNA pool by any RNase producing 2',3'-cyclic phos- 
phate, such as Irel, leaving only the unmodified pool as 
a possible substrate. It is tempting to speculate that the 
tRNA modification evolved to protect tRNA phe from 
cleavage by Irel and that regulation of the methylated 
pool of tRNA phe could provide a control point by which 
Irel may affect the translational capacity of cells. 

Conclusions 

We investigated the mechanism of RNA cleavage by 
Irel RNase by combining computational, structural and 
quantitative biochemical approaches. This work experi- 
mentally determined the location of the catalytic center 
in Irel RNase and converged on a model of docking 
and cleavage of cognate RNA stem-loops. Histidine 
H1061 was defined as a key catalytic residue that is 
required in single copy in the twofold symmetric active 
site of the Irel RNase. Presented here in trans mutagen- 
esis studies using (R1039A) and (R1039A, H1061N) Irel 
mutants, as well as electron density-guided in silico 
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docking of an RNA hairpin with the scissile bond cor- 
rectly aligned for catalysis suggest that the HLEs of two 
Irel monomers participate in stem-loop binding. This 
model intuitively predicts that only one RNA stem-loop 
binds to the dimer of Irel at a time (Figure 5e) whereas 
the high-order oligomerization of Irel enhances RNA 
recognition via allosteric stabilization of the HLE [5]. 
These findings identify the link between oligomerization 
of Irel and recognition of cognate RNA stem-loops, 
providing a glimpse of the central event in Irel 
signaling. 

Methods 

Experimental errors 

Individual rate constants were determined from single- 
exponential fitting of time courses of RNA cleavage. 
Kinetic parameters were reproduced two or more times 
and were consistent between different days. Rate varia- 
tions were typically within twofold; this uncertainty is 
small compared to the effects we describe as significant. 
Experimental errors are provided in the text and figure 
legends, when appropriate. 

Irel expression and purification 

IrelKR32 and its mutants were expressed as GST fusion 
proteins using pGEX-6P-2 plasmid (GE Healthcare, 
Waukesha, WI) and codon-compensated Escherichia coli 
(BL21-CodonPlus RIPL) (Stratagene, Santa Clara, CA). 
Expression was performed at room temperature for 4 h 
after IPTG induction. Cells were lysed using an Emulsi- 
flex C-3 (Avestin Inc. Ottawa, Ontario, Canada) and 
proteins purified by affinity purification with subsequent 
cleavage of the GST tag by Prescission protease (GE 
Healthcare, Waukesha, WI). All protein mutants were 
fractionated by gel filtration on an S200 column to 
approximately 99% purity. Proteins stocks (10-20 mg/ 
ml) were stored at -80°C in the presence of 5% glycerol. 

RNA substrates 

RNA oligonucleotides were purchased from Dharmacon 
(Lafayette, CO), labeled at the 5'-terminus using T4 
polynucleotide kinase and 32 P-ATP (PerkinElmer, Wal- 
tham MA) and purified by 20% polyacrylamide gel elec- 
trophoresis (PAGE) that allowed a single-nucleotide 
resolution, as described previously [5]. 

Irel RNase cleavage assay 

Kinetics of RNA cleavage was conducted as described 
previously [5], Typically, reactions were carried out in 
10 ul volume at 30°C. Reactions were started by adding 
1 ul of 32 P-labeled RNA to 9 ul of premixture contain- 
ing 20 mM 4-(2-hydroxyethyl)-l-piperazine-ethanesulfo- 
nic acid (HEPES) pH 7.4, 70 mM NaCl, 2 mM MgCl 2 , 4 
mM dithiothreitol (DTT) and 5% glycerol. The reactions 



contained < 1 pM radioactively ( P) labeled RNA and 
were conducted under single-turnover conditions 
(except for measurements of k cat and K m on Additional 
file 1, Figure SI, which were conducted under multiple- 
turnover conditions). Unless noted otherwise, Irel con- 
centration was 3 uM. The enzyme concentration was 
determined from IrelKR32 sequence using absorbance 
at 280 nm 

(e 28 o = 40.8 x 10 3 M/cm was calculated with 
BiochemLabSolutions ELN software, Princeton, NJ). 
Reactions were quenched at time intervals with 6 ul 
stop solution containing 10 M urea, 0.1% SDS, 0.1 mM 
ethylenediaminetetra-acetic acid (EDTA), 0.05% xylene 
cyanol and 0.05% bromophenol blue. Samples were ana- 
lyzed by 10% to 20% PAGE, gels were scanned using 
Typhoon (Molecular Dynamics-GE Healthcare, Wauke- 
sha, WI) and quantified using ImageQuant (Molecular 
Dynamics-GE Healthcare, Waukesha, WI) and Gel- 
Quant. NET (BiochemLabSolutions, Princeton, NJ) pro- 
grams. The data were plotted and fit in SigmaPlot to 
exponential curves to determine observed rate constants 
and to hyperbolic curves to determine binding 
constants. 

Irel oligomerization assay 

Irel was in the reaction buffer containing 20 mM 
HEPES (pH 7.4), 70 mM NaCl, 2 mM MgCl 2 , 2 mM 
ADP, 4 mM DTT and 5% glycerol. Transparency of the 
samples containing higher concentrations of Irel visibly 
changed immediately upon adding the enzyme. Samples 
were allowed to sit for 15 min to allow for a complete 
Irel oligomerization. The optical density of the samples 
was subsequently measured at room temperature (22°C), 
at 500 nm on a UV-visible spectrophotometer Ultrospec 
3300 Pro (Amersham-GE Healthcare, Waukesha, WI). 
OD 50 o was obtained after subtraction of baseline absor- 
bance from the buffer free of Irel. Sample absorbance 
did not change upon repeated scans of the same sample, 
indicating that the oligomerization reaches equilibrium 
under the conditions of the experiments. 

Irel-(dCdCdGdCdAdG) complex crystallization 

Irel-(dCdCdGdCdAdG) complex was prepared by mix- 
ing IrelKR32 and dCdCdGdCdAdG (IDT). When per- 
formed at relatively low NaCl concentration (300 mM 
or less) addition of the oligonucleotide causes profound 
precipitation of the oligomeric complex, as observed 
upon addition of ADP-Mg [5]. The oligonucleotide thus 
apparently promoted formation of IrelKR32 oligomers 
analogous to ADP-Mg and to APY29. The precipitate 
readily redissolved upon slight increase in NaCl concen- 
tration, indicating that the oligomerization reaction was 
salt dependent and readily reversible. Crystallization was 
conducted in hanging drops using stock solution of pre- 
made complex containing IrelKR32 (12 mg/ml) and 0.6 
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mM 2'-deoxy-CCGCAG. The well solution contained 

0. 12 M sodium citrate (pH 6.5), 7% polyethylene glycol 
(PEG) 3350 and 4% glucose. Single crystals grew over- 
night and were cryoprotected in well solution containing 
25% ethylene glycol. The crystals belong to orthorhom- 
bic space group C222 distinct from the orthorhombic 
space group previously reported for IrelKR32 
oligomer with bound APY29. Crystallization of a variety 
of other nucleotides and protein constructs was also 
tested but produced either inferior crystals or no 
crystals. 

Diffraction from Irel-(2'-deoxy-CCGCAG) crystals was 
collected on Beam Line BL 8.3.1 at the Advanced Light 
Source (Lawrence Berkeley National Laboratory, Berke- 
ley, CA), at an X-ray wavelength of 1.115872 A and an 
oscillation angle of 1°. The data were indexed, integrated 
and scaled using the XDS package [30] (Additional file 

1, Table SI, Figure S5). A total of 5% of the reflections 
were marked as a test set (Rf ree )- A molecular replace- 
ment solution was found using PHASER [31] starting 
from monomer C of PDB ID 3fbv as a search model. 
Seven copies of monomer C were found in the asym- 
metric unit and the resulting 7-mer of IrelKR32 was 
used for rigid body refinement in Phenix [32]. Simulated 
annealing was attempted and, expectedly, resulted in 
deterioration of statistics and excessive separation 
between R wor k and Rf ree (Additional file 1, Table SI). 
Simulated annealing was therefore used only for calcula- 
tion of omit maps (2000°C). Fourier c>A-weighted [33] 
Fobs-Peak difference maps were used for interpretation of 
the parts of the model missing from the starting struc- 
ture. Electron density and structure were analyzed and 
graphed in Coot [34] and PyMol (Schrodinger, San- 
Diego, CA). The sevenfold non-crystallographic symme- 
try (NCS) of the model was used to considerably 
increase quality of the electron density maps by aver- 
aging. However, all elements of Irel secondary structure 
were clearly visible without NCS. Coordinates have been 
deposited with the Protein Data Bank http://www.rcsb. 
org/ under accession number 3SDM. 

Crystallization of Irel KR32(H1 061 N) oligomer 

The H1061N mutant was crystallized and analyzed as 
described for the 3.2-A structure of Irel [5]. Diffraction 
data were collected on the beamline 8.3.1 at the 
Advanced Light Source (ALS, Berkeley, CA) using wave- 
length 1.115879 A. A total of 5% of reflections were 
marked as test set. The 3.2-A structure [5] was used for 
molecular replacement. The H1061N mutation was 
introduced in Coot [34] prior to refinement. Rigid-body 
refinement and subsequent simulated annealing refine- 
ment (starting temperature 3200 K) using 14-fold non- 
crystallographic symmetry were performed in Phenix 
[32] (Additional file 1, Table S3). The resulting crystal 



structure had a good stereochemistry and no residues in 
disallowed regions on the Ramachandran plot. In order 
to reduce model bias, we also calculated an omit differ- 
ence map from a model refined at 3200 K without the 
use of NCS restraints, from which a whole a-helix con- 
taining the H1061N residue was omitted prior to refine- 
ment. The resulting electron density clearly shows that 
the conformation of the H1061N mutant is identical to 
that of wild-type Irel. Coordinates have been deposited 
with the Protein Data Bank http://www.rcsb.org/ under 
accession number 3SDJ. 

Additional material 



Additional file 1: Gepasi modelling files. Compressed archive with 
Gepasi kinetics modeling files used in Figure S3. 

Additional file 2: Supplementary information. File with supplementary 
figures and tables. 
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